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Abstract 


A nanofluid is a suspension of ultrafine particles in a conventional base fluid which tremendously 
enhances the heat transfer characteristics of the original fluid. Furthermore, nanofluids are expected 
to be ideally suited in practical applications as their use incurs little or no penalty in pressure drop 
because the nanoparticles are ultrafine, therefore, appearing to behave more like a single-phase fluid 
than a solid—liquid mixture. About a decade ago, several published articles focused on measuring and 
determining the effective thermal conductivity of nanofluids, some also evaluated the effective 
viscosity. There are only a few published articles on deriving the forced convective heat transfer of 
nanofluids. The purpose of this article is to summarize the published subjects respect to the forced 
convective heat transfer of the nanofluids both of experimental and numerical investigation. 
© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Improvements to make heat transfer equipment more energy efficient would need to 
focus on miniaturization on the one hand and an astronomical increase in heat flux on the 
other. Heat transfer fluids such as water, mineral oil and ethylene glycol play a vital role in 
many industrial processes, including power generation, chemical processes, heating or 
cooling processes, and microelectronics. The poor heat transfer properties of these 
common fluids compared to most solids is a primary obstacle to the high compactness and 
effectiveness of heat exchangers. The essential initiative is to seek the solid particles having 
thermal conductivities several hundreds of times higher than those of conventional fluids. 
An innovative idea is to suspend ultrafine solid particles in the fluid for improving the 
thermal conductivity of a fluid. Many types of particle, such as metallic, non-metallic and 
polymeric, can be added into fluids to form slurries. However, the usual slurries, with 
suspended particles in the order of millimeters or even micrometers may cause some severe 
problems. The abrasive action of the particles causes the clogging of flow channels, erosion 
of pipelines and their momentum transfers into an increase in pressure drop in practical 
applications. Furthermore, they often suffer from instability and rheological problems. In 
particular, the particles tend to settle rapidly. Thus, although the slurries give better 
thermal conductivities, they are not practical. 

The use of particles of nanometer dimension was first continuously studied by a research 
group at the Argonne National Laboratory around a decade ago. Choi [1] was probably 
the first one who called the fluids with particles of nanometer dimensions ‘nano-fluids’. 
Compared with suspended particles of millimeter-or-micrometer dimensions, nanofluids 
show better stability and rheological properties, dramatically higher thermal conductiv- 
ities, and no penalty in pressure drop. Several published literature have mainly focused on 
the prediction and measurement techniques in order to evaluate the thermal conductivity 
of nanofluid. Recently, Trisaksri and Wongwises [2] reviewed the literature on the general 
heat transfer characteristics of nanofluids. It was found that only a few papers have 
discussed on the convective heat transfer of nanofluids. The purpose of this article is to 
review the literature mentioning convective heat transfer of nanofluids, including the 
experimental and theoretical investigations. 


2. Experimental approach 


A several years ago, the nanofluid has been found to be an attractive heat transport 
fluids. It has exhibited a significant potential for heat transfer augmentation relative to the 
conventional pure fluids. It has been expected to be suitable for the engineering application 
without severe problems in pipeline and with little or no penalty in pressure drop. 
However, the earlier researches revealed that most studies of the transport properties of the 
nanofluid were focused on the thermal conductivity. The convective heat transfer of 
nanofluids has received comparatively little attention in literature. There are few 
publications dealing with the convective heat transfer of nanofluids. The most productive 
studies have been continuously carried out by the following researchers. 
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cross section area, m? 

concentration of liquid solution, kg,/kg..) 
specific heat, J/kg K 

nanoparitcle diameter, m 

inner diameter of the tube, m 

acceleration of gravity, m/s? 

heat transfer coefficient, W/m? K 

humidity ratio of the air, kgy/kg, 

thermal conductivity of the nanofluid, W/m K 
length of the test tube, m 

mass flow rate, kg/s 

Nusselt number 

pressure drop, Pa 

Peclet number 

Prandtl number 

heat flux, W/m? 

heat transfer rate, W 

radius of tube, m 

radius of disc, m 

Reynolds number 

perimeter, m 

thickness, m or time in the Lattice Boltzmann method 
temperature, °C 

mean velocity, m/s 

overall heat transfer coefficient, W/m? K 
velocity in y-direction, m/s 

volume, m? 

velocity in z-direction, m/s 

distance along the tube, m or x-coordinate 
y-coordinate 

z-coordinate 

thermal diffusivity, m/s 

thermal dispersion coefficient, N/m* K 

rate of dissipation (for k—-e turbulent model) 
friction factor 

tangential coordinate, rad 

volume fraction 

inclination angle, degree 

thermal diffusivity coefficient 

mass diffusivity, m?/s 

turbulent kinetic energy (for k—e turbulent model) 
density, kg/m? 

tensor of the dispersed thermal conductivity 
viscosity, m/s 
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Subscript 

a annulus or air 
d particle 

des desiccant 

D dispersion effect 
f fluid 

hf heating fluid 

i inside 

in inlet 

Im log mean 

m mean 

(0) out 

out outlet 

nf nanofluid 

r radial direction 
T turbulent 

w tube wall 

x axial direction 


Li and Xuan [3], and Xuan and Li [4] presented an experimental system to investigate 
the convective heat transfer coefficient and friction factor of nanofluids for laminar and 
turbulent flow in a tube. The sample nanofluid used in this study was deionized water with 
a dispersion of Cu particles with below 100 nm diameter. The nanofluids, with different 
particle volume percentages, were used to determine the effect of the nanoparticle 
concentration on the heat transfer coefficient, the percentages being 0.3, 0.5, 0.8, 1, 1.5 and 
2%. The Reynolds number of the nanofluids varied within the range of 800-25,000. The 
experimental apparatus was mainly composed of a reservoir tank, a pump, a pipeline, a 
test section, a cooler, and a fluid collection tank. Test section was a straight brass tube with 
10mm inner diameter and 800 mm length and a 3.5 kW electric heater was used to obtain a 
constant wall heat flux boundary condition. To prevent the aggregation of the 
nanoparticles in the base fluid, a fatty acid was used as a dispersant to cover the 
nanoparticles. The experimental apparatus of their work was depicted in Fig. 1. 

The heat transfer coefficient of the nanofluids is calculated from the following equations 


hat D 
Nurt = 1 
ee (1) 
q 
jy = > 2 
Inf Ty — T; (2) 


The experimental results illustrated that the convective heat transfer coefficient of the 
nanofluids varied with the flow velocity and volume fraction as well as being higher than 
the base fluid at the same conditions. Compared with water, the Nusselt number of the 
nanofluid with 2% volume fraction of Cu particles was 60% higher. Furthermore, 
the Dittus-Boelter equation was not valid for the prediction of the Nusselt number of the 
nanofluids at various volume fractions. From the experimental data, new heat transfer 
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Fig. 1. The experimental apparatus of Xuan and Li [3] and, Li and Xuan [4]. 


correlations for the prediction of the heat transfer coefficient of nanofluids flowing in a 
tube could be established in the following form: 
For laminar flow 


Nuns = 0.4328 (1.0 + 11.2856% Pe}!8) Red§ Prof (3) 


nf 


For turbulent flow 


Nuns = 0.0059 (1.0 + 7.6286" 6586 Pe900!) Re? 2238 ppd (4) 


where Nu is the Nusselt number, Pe is the Peclet number, Re is the Reynolds number, and 
Pr is the Prandtl number. 


d 
Peg = (5) 
Anf 
D 
Repr = (6) 
Unf 
and 
Pee (7) 
One 


In order to calculate these parameters, the thermal diffusivity of the nanofluid was 
defined as 


knf — knf 
(PCp) ap d- )(PCp) ¢ gy b(PCp) 


In the case of the pressure drop, the results indicate that the friction factor of the 
nanofluid was almost equal to those of water under the same flow velocity and did not vary 


Anf = 


(8) 


802 W. Daungthongsuk, S. Wongwises / Renewable and Sustainable Energy Reviews 11 (2007) 797-817 


with the volume fraction. This implies that the nanofluids cause no penalty to pump power 
and the single-phase friction factor correlation can be used for more accurate predictions. 
The friction factor of the nanofluids was determined from the following equation: 


= APD 2g 


j 
ni L 


(9) 


Yang et al. [5] reported experimental results which illustrated the convective heat 
transfer coefficient of graphite nanoparticles dispersed in liquid for laminar flow in a 
horizontal tube heat exchanger. Nanoparticles of different sources were used in this study 
and focused on an aspect ratio (l/d) of about 0.02 (like a disc), because the addition of large 
aspect ratio particles into a fluid may dramatically increase viscosity compared with the 
single-phase. In this study, two series of nanofluids with different base fluids were used. 
The experiments were performed within the following ranges: the volume flow rates were 
between 62-507 cm?/min, the Reynolds number varied between 5-110, and the heating 
fluid temperature between 50-70°C. This article intended to study the effects of the 
Reynolds number, volume fraction, temperature, nanoparticle source, and type of base 
fluid on the convective heat transfer coefficient. The experimental system used in this study 
was shown schematically in Fig. 2. 

The following calculation procedure was used to evaluate the heat transfer coefficient of 
the nanofluids. 

The heat transfer rate into the nanofluids is 


Qat = MCp(Tout — Tin) (10) 
The heat transfer from the heating fluids is 
One = Moe CPne(Tou — Tindar (11) 


The energy balance between the heating fluid and the process fluid were differences 
around 25%. This is common for a small-scale testing equipment. 
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Fig. 2. Schematic diagram of the experimental apparatus of Yang et al. [5]. 
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The overall heat transfer resistance relationship is used to determine the heat transfer 
coefficient of the nanofluids, Anf: 
1 1 Də, Do 1 


= l 12 
U ha (A/A) 2k "Di he (12) 


The overall heat transfer coefficient can then be determined from 
L 
AAT im 


The Monrad and Pelton’s equation for turbulent flow in annuli was used to calculate the 
outside heat transfer coefficient, Ao, as follows: 


U= (13) 


0.53 


Nu = 0.020Re°*Pr'/3 (D2/D1) (14) 


Subscript 2 and 1 mean the outer and inner diameter of the annulus, respectively. The 
equivalent diameter, De, was used instead of the typical diameter, D. 
2 2 
D.= Di, a Ds 
e D, 


The experimental results illustrated that the heat transfer coefficient increased with the 
Reynolds number and the particle volume fraction. For example, at 2wt% the heat 
transfer coefficient of the nanofluids was moderately increased compared with the base 
fluid, but at 2.5 wt% it was 22% higher than the base fluid at 50 °C fluid temperature and 
15% higher at 70°C. Several factors may lead to smaller improvements in heat transfer 
enhancement of nanofluids with increased fluid temperature. One is the rapid alignment of 
nanoparticles in less viscous fluids, leading to less contact between the nanoparticles. A 
second factor is that depletion of particles in the near-wall fluid phase could lead to a lower 
thermal conductivity layer near the wall. Two graphite nanoparticle sources at same 
particle loading give different heat transfer coefficient, as in the case of different base 
fluids. This may have been caused by particle shape, morphology or surface treatment. 
This means that one type of nanoparticle is more effective than another in augmenting the 
heat transfer coefficient of the nanofluid. Moreover, the correlation established by Li and 
Xuan [4] for laminar flow of nanofluids gave a higher heat transfer coefficient than that 
calculated from the experimental data. This may be because the effect of particle aspect 
ratio was not included. Similarly, the single-phase correlation of laminar flow in tubes 
failed to predict the convective heat transfer coefficient of nanofluids. 

All of the experimental data were used to develop a new heat transfer correlation for the 
prediction the heat transfer coefficient of laminar flow nanofluids in a more convenient 
form by modifying the Seider-Tate equation [6] as follow: 


aie pA T 
Q = Nun Pra” (5) (=) = aR, (16) 
wW 
where a and b are the various constants which depend on the experimental conditions. 
Wen and Ding [7] presented an experimental study which evaluated the convective heat 
transfer coefficient of yAl,O3 nanoparticles suspended in deionized water for laminar flow 
in a copper tube, focusing in particular on the entrance region. The various concentrations 
of nanoparticles (0.6, 1.0, and 1.6%) were tested under a constant wall heat flux. The 
experimental apparatus mainly consisted of a straight copper tube with 970mm length, 


(15) 
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4.5+0.02 mm inner diameter, and 6.4+0.05mm outer diameter as a test section. A 300 W 
silicon rubber flexible heater was used to keep a constant heat flux boundary condition. A 
peristaltic pump was used to deliver a maximum flow rate of 101/min. The flow rate was 
controlled by adjusting rotational speed of the pump. In order to stabilize the 
nanoparticles, sodium dodecylbenzene sulfonate (SDBS) was used as a dispersant. Their 
experimental apparatus was shown schematically in Fig. 3. 

The heat transfer coefficient and the Nusselt number of the nanofluids were calculated 
from the following equations: 


lop xD 
Nuns = (17) 
i kof 
q 
h = 18 
Inf,x Tux = Tex ( ) 


The energy balance was used to determine the fluid temperature profile along the test 

tube as follows: 
qSx 

Trx = Tin + — 19 

f,x wt Cad ( ) 

The results showed that the local heat transfer coefficient varied with the Reynolds 

number and volume fraction. In particular, the entrance region illustrated that the use of 

nanofluids affected a pronounced increase in the heat transfer coefficient causing a 

decrease in the thermal boundary layer thickness and will decrease along the tube length. 
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Fig. 3. Schematic diagram of the experimental apparatus of Wen and Ding [7]. 
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This behavior implied that it might be possible to create a ‘smart entrance’ region to meet 
the highest performance of nanofluids. Furthermore, the Shah correlation for laminar flow 
and the Dittus-Boelter equation for turbulent flow do not coincide with the Nusselt 
number between the experimental results and the calculated values. 


3. Mathematical approach 


Similar to the experimental investigation, many researchers devoted to formulate the 
appropriate correlations for prediction of the thermal conductivity of the nanofluid. It is 
difficult to establish any formulated theory for the prediction of the accurate heat transfer 
behavior of the nanofluid because the several complex phenomena may coexist in the main 
flow e.g. sedimentation, Brownian force, dispersion, and friction force between the fluid 
and the solid particles. Some researchers treat the nanofluid more like a fluid rather than a 
liquid—solid mixture. There are few articles taking the nanofluid into account as the multi- 
phase feature. Some detailed reviews mentioned above are shown as follows: 

Xuan and Li [8] presented a method for preparing some nanofluid samples and 
introduced a theoretical study of thermal conductivity. They used the transient hot-wire 
method for measurements. They then suggested a dispersion model or single-phase model 
to describe the heat transfer performance of nanofluids flowing in a tube. The 
nanoparticles used in this article were Cu particles with 100nm diameter, respectively, 
suspended in water and mineral oil. Oleic acid and laurate salt were used as dispersants in 
order to stabilize the particles in the mineral oil and water, respectively. For Cu—mineral 
oil, the results showed that the suspension reached a stabilized state within about 1 week 
and had no sedimentation with a fill of 22 wt%. In the case of Cu—water nanofluid 
suspension with 9 wt% laurate salt, vibrated by a ultrasonic vibrator, the results indicated 
that the nanoparticles can suspend more 30h than in stationary state. This behavior 
implied that the fluid viscosity might be an important factor affecting the dispersion and 
stability of the particles. Moreover, this study illustrated that the Hamilton and Crosser 
model [9] and Wasp model [10] for predicting thermal conductivity can be applied to 
obtain an approximate solution in the absence of a sophisticated formula. Using the 
transient hotwire method to measure the thermal conductivity, the results showed that the 
thermal conductivity remarkably varied with the volume fraction of the particles. There 
are two different approaches to analyze the heat transfer enhancement. One is the two- 
phase model and other is the single-phase model. Although the first model describes an 
understanding of the functions of both liquid phase and solid phase in the heat transfer 
process, it needs a long time for computation and a high performance computer. The 
second model, assuming that both the liquid and particle phases are in thermal equilibrium 
and flow at the same velocity, is simpler and needs less computation time. The purpose of 
their work aimed to develop a modified single-phase model to describe the heat transfer 
process of nanofluids flowing in tubes. The governing equations were expressed as follow: 


oT | OT | bB Sor. 19 b. ôT 
at ox = (0 = oa) ax FOr (0 i a a os 


If the axial temperature distribution is ignored, Eq. (20) can be reduced as 


ôT aT 12a By oT 
at’ “Ox For (Grk ra 
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The boundary conditions for this situation are as follows: 
T\,2r = Ty and T|,~9 = To (22) 


The application of the separation of variables method for plug flow, the following 
solution was obtained: 


—02,%/Pe JolOm?) 
=e ye J (Om)Om (3) 


where F= r/R, X= x/L, Nu=2Rh/k*,, Pœ = uL/ġ*, Pe = Pe* (R/L) and eigenvalues 
Om s are the positive roots of the following equation: 


T-T, 


JOm) = 0 (24) 
The heat balance can be expresses as 
T 
se Q =h(T»y — Ty) (25) 
or 1=R 


The analytical correlation of the Nusselt number is then given by 


oo eo 2 X/Pe 
nmi? (26) 
Saar Om / Pe J02, 


m= 


Nu= 

For presence of the dispersion in axial and radial direction, Eq. (20) and the modified 

Dankwerts’s boundary conditions (Eqs. (27a) and (27b)) in the axial direction and the 
boundary condition in radial direction Eq. (22) give solution (Eq. (28)). 


jet 


nf ay T = uApCp(T) — T) at x = 0, (27a) 
OF =o Geet (27b) 
Ox 
T= Tw `. X(x) J (On?) 
— "~ -2 28 
T T, 2 KO- LOPE] On a 


where X(X) and X’(X) are defined as follows: 


X (x) = ment S me” +™2* 


X'(%) = mm (etm = erie) 


Pet + y/ Pe? + 402, (L/R) 
2 
The eigenvalues are described in Eq. (24). Thus, the Nusselt number correlation for the 
nanofluids is obtained from 
°°, X(x)/|X(0) — X’(0)/Pe* 
Nu = LXO XO -XO Pe] m 
Dra LX(%)/[X(0) — XO /Pe]] /O;, 


m2 = 
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It should be noted that Eqs. (26) and (29) are derived for constant wall temperature and 
need the iteration procedure to solve the problem. The Nusselt number for the other 
boundary conditions can derive with the same method. 

Finally, this study expressed that either theoretical expressions or experimental work 
needed to explain the heat transfer process of nanofluids, in which several phenomena may 
coexist in the main flow. 

Xuan and Roetzel [11] presented two different approaches to determine some 
fundamentals to predict the convective heat transfer coefficient of nanofluids under 
the assumption that they behaved like a single-phase liquid rather than a normal 
solid—liquid mixture. They included both the effects of transport properties and thermal 
dispersion of the nanofluid. The first model treated the nanofluid as a single-phase 
fluid and the second model as a multi-phase and dispersed fluid. In a conventional 
approach, the assumptions were no slip velocity between the particles and liquid phase 
and that they were in a thermal equilibrium state. The nanofluid behaved like a common 
pure liquid. This implied that the whole equation of continuity, motion and energy for a 
single-phase fluid had been applied directly to the nanofluid. Similar to the Xuan and Li 
model [8], the Hamilton and Crosser model [9] and Wasp model [10] for predicting the 
thermal conductivity of nanofluid give rough estimations. For a modified conventional 
approach, called the dispersion model, the slip velocity between the solid—liquid mixture 
may not be zero causing many factors, such as Brownian force, gravity, friction force 
between the solid—liquid mixture, Brownian diffusion, sedimentation, and dispersion to 
happen simultaneously in the main flow. The details of two approaches are expressed 
as follows: 

For the conventional approach, the general form of the energy equation for the 
incompressible pure fluid with negligible viscous dissipation is 


T 
a+ V-uT =V-(aVT) (30) 


This equation is suitable to explain the heat transfer process of the nanofluid. This 
means that the heat transfer equation for the single-phase fluids can extend to the 
nanofluid directly. It should be noted that the thermal properties used in this equation are 
those of nanofluids. For this reason, the dimensionless heat transfer correlations of the 
single-phase fluid can be applicable for the nanofluid, so 


Nu = 3.66 (For fully developed laminar flow in tube, Tw = const.) (31a) 


Nu = 0.023Re°*Pr'/? (For turbulent flow in tube) (31b) 


For the modified conventional approach, the irregular movement of the ultrafine 
particles led to little disturbances of both the temperature, 7’, and velocity of the 
nanofluid, u’. Hence, the real temperature and velocity averages are expressed as 


T=(T¥ +T (32a) 


u= (uy +u (32b) 
where (TY = (1/V¢) Sy, TAV, (uy = (1/ V9 Sy, wdV and (1/V1) fy, T'dV = 0 
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By the procedure demonstrated by Kaviany [12] and assuming that the boundary surface 
between the fluid and the particles is so small, Eq. (30) can be rewritten as 


OT)’ 
(PCy) ap ôt 


+ (uY V- (TY | = V- (kar YTY) — (PCp) pV T'Y (33) 


The chaotic movement of the ultrafine particles in the flow gives the thermal dispersion 
(the second term on the right hand side of Eq. (33)). Therefore, the heat flux resulted from 
the dispersion in the nanofluid flow can be expressed as: 

dp = (PCP) (u T'Y = =p - V(TY (34) 
where t is the tensor of the dispersed thermal conductivity. Thus, Eq. (33) can be 
rewritten as 

nae totum | =¥- e VY 
+ (uY V- (TY | =V. || anI + V(T'Y 35 
ôt ( ) ( ) nf oT ( ) ( ) 


Then, the energy balance equation for the nanofluid can be given as 


/ 


OT 
(OCP) + PCat) YIT) = (PCP) YO ITY) G6) 


where y may be called the total effective thermal diffusivity tensor. It combines the 
molecular effect and the thermal dispersion effect. y is defined as follows: 


knf eT 
(p Cp Jaf 


It is anticipated that the thermal dispersion tensor is a function of the flow pattern, 
properties of the base fluid and the nanoparticles, dimensions and shape of the 
nanoparticles, and volume fraction of the nanoparticles. 

The energy equation for the nanofluid flowing inside a tube for this situation may be 
expressed as 


y= 


+(1— Wp (37) 


oT oT 10 kp, oT 
ôt E ôx rér | (aa % ze) Z i 
For taking the thermal dispersion in the axial direction into account, so 
ôT oT 10 kp, oT ô kpx \ ôT 
= n : nf i FT 39 
at ax ror (e = a d Tax (e a ae d ea 


The thermal diffusivity coefficients can be found from the derivation of Taylor [13] and 
Aris [14]. Unfortunately, there is no theoretical and experimental work established on the 
dispersed thermal diffusivity of the nanofluid. The kp appearing in Eqs. (38) and (39) can 
be calculated by Plumb [15] and Hunt and Tien [16]. 

To take the effect of thermal dispersion into account, the surface heat flux can be 
expressed as 


q=- (har Sot hor Sean (40) 


Similar to the procedure proposed by Nield and Bejan [17], the dimensional analysis of 
the above equation formally resulted in the heat transfer correlation of the nanofluid as 


W. Daungthongsuk, S. Wongwises / Renewable and Sustainable Energy Reviews 11 (2007) 797-817 $809 


follows: 

Nux = [1 + C* Pef’(0)]0'(0)Re” (41a) 
with considering the flow pattern, Eq. (41a) can be modified to Eq. (41b) 

Nux = [1 + C*Pe"f’(0)]0'(0)Re” (41b) 


In which Pe = (rujy/Onp), C’ is an unknown constant determined from experiment, m and 
n depend on the flow pattern, fis dimensionless velocity, 0 is dimensionless temperature, 
and f’ and 6’ are the derivatives of the dimensionless velocity and dimensionless 
temperature, respectively. 

Maiga et al. [18] presented a mathematical formulation and numerical method to 
determine the forced convective heat transfer and wall shear stress for the laminar and 
turbulent regions of water—yAl,O3 and ethylene glycol-yAl,O3 flowing inside a uniformly 
heated tube. The solid—liquid mixture takes a single-phase behavior into account, so the 
slip velocity between the phases was neglected. Furthermore, the local thermal equilibrium 
of the mixture and symmetry in flow were considered. For the turbulent flow region, the 
Reynolds-averaged Navier-Stokes equation and k—e turbulent model were used to describe 
the stresses and heat flux of the nanofluids. The conditions of the numerical simulations 
were as follow: the heated tube had 0.01 m diameter and 1.0m length. In the situation of 
laminar flow, the Reynolds number was fixed at 250 and the constant heat flux was varied 
between 10-250 W/m”. For turbulent flow, the constant heat flux was fixed at 500,000 W/ 
m? and the Reynolds number was varied in the range of 10,000-50,000. The problem 
configuration for this study is shown in Fig. 4. 

Under the above conditions the governing equations depicted in a (R, œ and Z) 
coordinate system are as follows: 


V- (pu) = 0 (42) 
V - (puu) = —OP/0X; + V - (uVuj) + Si (43) 
V: (puCpT) = V - (KYT) (44) 


where u = (ur, Uo, Uz) is the velocity vector, X; represents a spatial direction (X; = R, œ 
and Z), and S; represents the remaining viscous term corresponding to the coordinate 
system. 


— for the radial direction, R (i= 1); 


PUugug uR 2 Ou 
Ss; = 
i= p ” E R = 


(45) 


Fig. 4. The problem configuration in the study of Maiga et al. [18]. 
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— for the tangential direction, œ (i = 2); 


2 TO TT 
S=] OuR a puru 


46 
R ôw R R (46) 
— for the axial direction, Z (i = 3); 
53 =0 (47) 


For the turbulent flow regime, this article shows the application of the famous semi- 
empirical k—e turbulent model. The equations for this case are described as follow: 


Vou) =V: [n+ ft vn] + 0x- pe (48) 
Prox 


Ur 
V. o) =V. 
(pue) e F Pres 


ve! + Ci(e/K) Gx, + Cap? /K (49) 


where G, is the generation of turbulence kinetic energy due to the mean velocity gradients, 
and ur is the turbulent viscosity calculated from 


Hr = pCa? /e (50) 
Various constants in the above equations are as follows: 
Cis = 1.44; Co, = 1.92; Cy, = 0.09; Prr = 1.0; and Prr = 1.3 (51) 


The numerical results illustrated that the heat transfer and the wall friction of nanofluids 
increased with increasing the particle concentration and the ethylene glycol-yAl,03 gave 
higher heat transfer enhancement than the water—yAl,O3. Then, the heat transfer 
performance of the nanofluids was more pronounced with the increase of the Reynolds 
number for the turbulent flow region. 

Roy et al. [19] proposed a numerical investigation to evaluate the heat transfer and wall 
shear stress for radial laminar flow in a cooling system of water—yAl,O3 nanofluids 
compared with water, glycol, and oil as base fluids. All assumptions used in this study were 
similar to Maiga et al. [18] who took the nanofluid as an incompressible fluid. The 
geometric configuration of this work is shown in Fig. 5. 


Z 
A 
= inlet SS 
ee 
iy 
Outlet <— <x la 
A 
Ro! 
| Symmetry of 
pe R | Revolution 


Fig. 5. The geometric configuration in the study of Roy et al. [19]. 
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The governing equations in the cylindrical coordinate (R, Z) and symmetry systems can 
be given as follows: 


V. (pu) =90 (52) 
V - (puu) = —VP + V-(uVuj)+ Si, i= 1,2 (53) 
V-(puCpT) = V - (kYT) (54) 


— for the radial direction, R(i = 1); 


Sau H (55) 


— for the axial direction, Z(i = 1); 


S2 =0 (56) 


The results show that the heat transfer rate and wall shear stress were higher with 
increasing particle concentration and Reynolds number. 

Nguyen et al. [20] introduced a numerical simulation to determine the efficiency of 
water—yAl,03 and ethylene glycol-yAl,O3 nanofluids for the cooling of a high-heat output 
microprocessor under laminar forced flow inside a heat sink. The simulation set-up 
consisted of a 50 x 50 x 10 mm rectangular slot with a 3 x 48 mm fluid flow cross-section. 
Assumptions used in this article were laminar flow, uniform velocity and temperature 
profile at the inlet section. The contact area for heat exchange was 10 x 10mm. The 
problem configuration for this study was sketched in Fig. 6. 

The numerical results indicate that the average heat transfer coefficient of nanofluids 
varied higher with the volume fraction and Reynolds number and the ethylene 
glycol-yAl,0;3 nanofluid gave higher heat transfer coefficient than the water—yAl,03 
nanofluid. 

Ali et al. [21] reported the mathematical and numerical formulation to determine heat 
and mass transfer between air and falling solution film in a cross-flow heat exchanger for 
the dehumidification and cooling process. The Cu nanoparticles were added into the 
solution in order to evaluate their effect in the augmentation of the heat and mass transfer 
process. Two models, proposed by Xuan and Roetzel [11], were used in this study. 
Furthermore, the assumptions of this article were as follow: (1) laminar and steady state 
flow; (2) thermal properties of the air and the solution are constant except for the thermal 


Contact Area with Outlet 
Microprocessor (10x 10 mm) 


aX <Q (Ww) ink 
y AA Heat sink (Copper Body) 


Through Flow Inlet Free Flow 
Cross-section 


Fig. 6. The problem configuration in the study of Nguyen et al. [20]. 
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Flow of Liquid 
Film Desiccant 
with Cu nanoparticles 


4 WANS 
tvyvvty 


Ydes <—+ 


Flow of Air 


OE 


A Te Z 7 8 >Ja 


Fig. 7. Schematic diagram of problem configuration of Ali et al. [21]. 


conductivity of the solution; (3) gravitational force of the air is neglected; (4) constant film 
thickness and (5) fully developed velocity profile for fluids. The schematic diagram of this 
work was shown in Fig. 7. 

The governing conservation equations of mass, momentum, energy, and mass diffusion 
are illustrated in two portions; air and solution (liquid film) as follows: 


— For the air 
Ow, 
a7 ~ 0 (57) 
OP, Ow, 
az = Hal aye ) (58) 
oT. OT, 
PaCPaWa “ez = ka ay2 (59) 
OH, @ Ca 
Va $ — Y 
W zz oy (60) 
— For the solution (liquid film) 
OUdes = 
= 0 (61) 
OUdes 
Paes + Hdes (=) =0 (62) 
Vies 
OT aes ð OT des 
Cc = —— | (kanf + Kais 63 
o Pues ox OY des [ os a C ) 
OC des om Cues 
a a. = es, 4 
ae ae (64) 
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The energy balance equation at interface between air and film desiccant is 


OT, OH, OT des 
ka — wih = —k 65 
a dy, + pala fg dy, des OVaes ( ) 
at Ya = ta, Ydes = taes 0<x <H and 0<z< L 

The mass balance at the interface is 


ôH a ð Cds 
= ( Y) es (66) 
oy a oy des 


a 
at Ya = ta, Ydes = taes OS X<HA and 0<z< L 
With the suitable boundary conditions, the velocities for both air and film desiccant are 
obtained analytically from Eqs. (58) and (62) as follows: 


1 OP, 


Waa) = 2, az 


P a Pa 


04-4) (67) 


tides aes) = EEST Yes toes — 7S) (68) 
Hdes 


The air side pressure drop and the solution thickness are calculated from 
OP, 3. Uaa 


éz 2Hp,B (oa) 
f 1/3 

TAE ts] (70) 
Pdes9 


The results illustrated that the dehumidification and cooling process increased with the 
following parameters: low air Reynolds number, high Cu particle volume fraction, 
increases in the height and length of the channel, and decreases in the channel width. Also, 
it was found that an increase in Cu volume fraction led to more stability of the solution. 

Ali and Vafai [22] presented the mathematical and numerical formulations to evaluate 
the effects of the inclination angel of parallel-and-counter-flow on heat and mass transfer 
between air and falling desiccant film, with Cu nanoparticles suspended. As Ali et al. [21], 
the same assumptions were held and two approach models, introduced by Xuan and 
Roetzel [11] where used in this literature. Then, the simulated conditions used in this study 
are summarized as follow: low and high air Reynolds number for the inclined parallel flow 
arrangement, and low and high air Reynolds number for the inclined counter flow 
configuration. The problem configuration for this study was shown schematically in Fig. 8. 

Similar to Ali et al. [21], the general governing equations (continuity, momentum, 
energy, and mass diffusion for air and the desiccant film) for both parallel flow and counter 
flow can be given by 


— For the air 
OUa  OVa 

eo 1 
ae Oi (71) 


Ôa Ôa oP ua 
a a = a 72 
(e etage) atalar) ms 
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Film Desiccant with Film Desiccant with 
Ultrafine Particles | Ultrafine Particles 


t, i 
eee y2 
Air Flow vx 


Air Flow y 


H 
(a) parallel flow (b) counter flow 
Fig. 8. The geometric configurations in the study of Ali and Vafai [22]. 
OV, OV, OP Ov, 
= 73 
ra(m set os) = atalas as 
OT... OT, oT, 
ne) a — | = K,| — 74 
Pa p(n ax, + Va =) ( D (74) 
ôH ôH OH 
= Y| — 75 
“aa t By, (a) ae 


— For the desiccant film 


OUdes 
=0 76 
Ox) ( ) 
@ udes 
pasg + Ho) =0 (77) 
OT aes 6) OT des 

= — |(k k is 

(PC rhs = ge (C + ha) S| (18) 
2; 

TR a = (79) 


Ria een pee 
ax» des dy 


In order to solve the governing equations of the air, the coordinate transformations are 
needed for both inclined parallel and counter flow. The simulation conditions used in this 
work are in the range of low and high air Reynolds number. 

From the numerical results it could be concluded that the dehumidification and 
cooling process of air and regeneration process of the liquid desiccant are enhanced with 
increasing inclination angle for both parallel-and-counter flow-channels. Moreover, the 
effective thermal conductivity of the liquid desiccant varied with the Cu nanoparticles and 
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the thermal dispersion and low desiccant concentration led to enhancement of the 
dehumidification process. 

Xuan et al. [23] proposed the thermal Lattice Boltzmann model to display the flow 
features and heat transfer process of Cu-water nanofluid flowing inside a channel. The 
important advantage of this method was that it took the molecular dynamics into account 
and bridged the gap of microscopic or macroscopic phenomena of the nanofluids. For 
application of this method to model the thermal fluid flow problem, there are two existing 
methods: the multispeed (MS) and the double-distribution-function (DDF), the latter 
method being used in this article. This led to two separate distribution functions as follow: 
the temperature is treated as a passive diffusing scalar, and it is simulated by a distribution 
function independent of density distribution, respectively. The feature point of the Lattice 
Boltzmann method is the assumption that the particles are mesoscopically located at a 
series of lattices and their distributions correspond to the Boltzmann distribution. For 
simulation of heat transfer characteristics of the nanofluid, the viscosity dissipation term in 
the Lattice Boltzmann equation for energy transport will be neglected. The equation is 
given by 

| Ti(x, ù — T(x, | | 


T(x + eAt,t + At) — T(x, t) = - > i=1,2,3,4 (80) 
F 


where e; is the lattice velocity vector, tr is the collision-relaxation time constant for the 
temperature distribution. T;" is the equilibrium function for the temperature which is 
defined as 


i T 
T= AT+ BEI, i=1,2,3,4 (81) 
where c is the reference lattice velocity (Ax/Ad) 

By using the Lattice Boltzmann method to model the nanofluid, the dynamic factors 
influencing the flow behavior of the nanofluid must be considered first. Without external 
forces, the nanoparticles are under influence of the buoyancy force, gravitational force, 
Brownian force, drag force and dispersion force resulting from the repulsive potential [13]. 
Hence, the sum of the total forces vector acting on the nanoparticle per unit lattice volume 
(V) is given by 

F=n(Fat+Fet+Fp+Fu)/V (82) 


where n is the number of the nanoparticles in a given lattice, F4 is the dispersion force, Fg 
is the Brownian force, Fp is the drag force, and Fy represents the buoyancy and 
gravitational force, respectively. These several variables are defined as 


8 
OV. 
Fas) on (83) 
=i fi 
1 2a? 2a? L — 4a? 
Fia Alt tn z ) (84) 


Fg = Giy C/At (85) 


C=2kgT; 4= 6ra (86) 
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Fp = —6zpaAu (87) 


4 
Fy = -370 gAp (88) 


where a is the radius of the nanoparticle, Ap and Au are the density and the velocity 
difference between the nanoparticle and the base fluid, respectively, 2 is the friction 
coefficient, kg is the Boltzmann constant, G; is the zero mean (unit variance independent 
Gaussian random number), L is the distance between center-to-center of the particles, and 
A represents Hamaker constant which accounts for the material properties but not for the 
configuration shape. 

The numerical results showed that, for 1% volume fraction the Nusselt number of 
Cu-water nanofluid was 27% higher than that of pure water at the same Reynolds 
number. Moreover, fluctuations of the Nusselt number of the nanofluid taking place along 
the main flow direction may cause random motion of the particles suspended in the fluid. 


4. Conclusion 


Evidently, nanofluids have great potential for heat transfer enhancement and are highly 
suited to application in practical heat transfer processes. This offers an opportunity for 
engineers to develop highly compact and effective heat transfer equipment. Several 
published articles show that the heat transfer coefficient of nanofluids are much higher 
than the common-base fluid and give little or no penalty in pressure drop. The main 
reasons for the heat transfer enhancement of the nanofluids may roughly be listed as 
follow: the suspended nanoparticles increase the thermal conductivity of the fluids, and the 
chaotic movement of ultrafine particles increases fluctuation and turbulence of the fluids 
which accelerates the energy exchange process. Furthermore, that the numerous 
correlations for single-phase fluid have clearly failed to predict the heat transfer 
coefficients of nanofluids may be caused by the complicated phenomena that coexist in 
the main flow. However, some sophisticated correlations for predicting the heat transfer 
performance of nanofluids have been established. At the present time, the uses of the 
nanofluids ar still remaining in early stages of development. Therefore, urgent theoretical 
and experimental research works is needed in order to clearly understand and accurately 
predict their hydrodynamic and thermal characteristics. Generally, many researchers 
indicated that nanofluids behave like pure fluids because the suspended particles are 
ultrafine. However, at present, no formulated advanced theory exists to explain the 
behaviour of nanofluids by considering them is as multi-component materials. 
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